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Abstract The standing slow magneto-acoustic oscillations in cooling coronal loops are investigated. 
There are two damping mechanisms which are considered to generate the standing acoustic modes 
in coronal magnetic loops namely thermal conduction and radiation. The background temperature is 
assumed to change temporally due to optically thin radiation. In particular, the background plasma is 
assumed to be radiatively cooling. The effects of cooling on longitudinal slow MHD modes is analytically 
evaluated by choosing a simple form of radiative function that ensures the temperature evolution of 
the background plasma due to radiation coincides with the observed cooling profile of coronal loops. The 
assumption of low-beta plasma leads to neglect the magnetic field perturbation and eventually reduces the 
MHD equations to a ID system modelling longitudinal MHD oscillations in a cooling coronal loop. The 
cooling is assumed to occur on a characteristic time scale much larger than the oscillation period that 
subsequently enables using the WKB theory to study the properties of standing wave. The governing 
equation describing the time-dependent amplitude of waves is obtained and solved analytically. The 
analytically derived solutions are numerically evaluated to give further insight into the evolution of the 
standing acoustic waves. We find that the plasma cooling gives rise to a decrease in the amplitude of 
oscillations. In spite of the reduction in damping rate caused by rising the cooling, the damping scenario 
of slow standing MHD waves strongly increases in hot coronal loops. 
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1. Introduction 


Observation made by SUMER instrument on board SoHO confirms that longitud i nal standing (slow 


magneto-ac oustic waves have been seen to suffer a rapid damping ( Kliem et ai . l200^ Wajig__gt__a4 


2002 200^ . Modelling for observed oscillations by SUMER have been carried out bv iTarovan et all 


(12007 1. where excitation and damping mechanisms are investigated. It is also presented how standing 
and propagating slow MHD waves can be differentiated from each other. The issue of fast damping is 
still obscure and attracts a remarkable attentions to understand the dominant method of damping. The 
most nominated mecha nism for damping the standing longitudinal oscillatio ns are thermal conduction, 
adiation and viscosity (Offnan and WansJ. 20021: De Moortel and Ho^. 20031 : Mendoza-Bricen o, Erdelyi, 


and Sigalottij2004j^ro^an_et^ 


2p05|;|Si galotti. Mendoza-Briceno, and Luna-Cardozo|, 200?]^ Al-Ghafri 


and Erdelyi, 201,01 : Al-Ghafri et oL . 2014ll . Radiation is found to dominate the damning of cool coronal 
loops whereas hot l oops are dam ped by thermal conduction. For a recent review of standing slow waves 
in coronal loops see I Wan Sdmil)- 


Plasma cooling has been detected everywhere in the solar atmosphere (jViall and Klimchuld . 12012 


201311. In the abs ence of coronal heating, the plasma starts to coo l by radiation and thermal conduction 
( KlimchuS . l2012[ l. For example. lAschwanden and Terradas ( 2008h have shown that the coronal loops are 
cooling with the characteristic cooling time of the order of a few oscillation periods. iMorton and Erdelyi 
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( 20091 . [ 201 ^ : Ruderman ( 2011a 3) investigated the effect of cooling on kink oscillations of coronal loop. 
They found that the kink oscillations experience an amplification due to the cooling. Further to this, 
propagatin g slow MHD waves in a homogeneous, radiativelv cooling plasma are studied bv Morton. Hood, 
and Erdelvi (|2nin ). Ess entially, plasma cooling is report e d to c ause a strong damping for longitudinal slow 
MHD waves. Recently. lErdelvi. Al-Ghafri. and MortonI ( 20111 1 investigated the behaviour of longitudinal 
magneto-acoustic oscillations in slowly varying coronal plasma. In particular, the damping rate is found 
tojmdergo_a^£edu£l^^ the emer gence of cooling. 

Al-Ghafri and Erdelyi ( 2013l l and lAl-Ghafri etal. ( 2014ll have studied the effect of cooling on standing 


slow MHD waves in hot coronal loops that are damped due to thermal conduction. In particular, the 
plasma cooling is assumed by unspecified thermodynamic source. The individual effect of cooling is found 
to amplify the amplitude of oscillating loops. Although the rate of damping caused by thermal conduction 
increases at the beginning, it is noticed to reduce gradually in a hot corona because of plasma cooling. 
As thermal conduction approaches infinity, the damping rate tends to zero. Hence, slow magnetosonic 
waves propagate almost without damping at the slower, isothermal sound speed. 

The present article discusses the combined effects of radiative cooling and thermal conduction on 
damping longitudinal standing MHD waves. A function with a simple form representing the radiation 
mechanism is taken to ensure that the plasma temperature decreases exponentially with time according 
to observations. The cooling is assumed to be weak with the characteristic cooling time scale much larger 
than the oscillation period. The paper is structured as follows. In Section d^]) we present our model and 
derive the main governing equation with boundary conditions. The analytical results are obtained with 
the aid of the WKB theory in Section ([31). In Section (|1]) the individual and combined effects of radiation 
and thermal conduction are studied by displaying the analytical solution numerically. Our discussions 
and conclusions are presented in Section d^. 


2. The Model and Governing Equations 


We model a straight coronal loop in which the magnetic field is uniform and in the z-direction, i.e. 
Bq = BqZ. The magnetic loop is of length L with its ends at z = ±L/2 . The background plasma 
temperature (pressure) is assumed to change as a function of time due to radiative cooling. The density 
is a constant. 

The governing MHD equations for the background plasma take the following form 


| + V,(pv) = 0, 

5v 1 

-I- p(v.V)v = -Vp -I -(V X B) X B, 

at po 


R 


M (7 - 1) 

P = —pT, 

M 

- = Vx(vxB). 


9 T , T 


= v(K||Vr)-p2Q(r), 


( 1 ) 

( 2 ) 

( 3 ) 

( 4 ) 

( 5 ) 


Here, v, B, p, p and T represent the flow velocity, magnetic field, density, pressure and temperature 
respectively; po is the magnetic permeability of free space, R is the gas constant, p is the mean molecular 
weight and 7 is the ratio of specific heats. The thermal conduction term is V(K||Vr) where K|| = 
and p^Q{T) is the general radiation term for optically thin losses. 


Note 

Tucker, 


that the radiativeloss,JunctionAs,^£groximate^^^^^^iecewise_continuous,^mctioi^Rosner, 
and Vaiana, ll97Sl: Ip ries tl. l2nn(ti but' it is assumed here to take a simple form to mach the 
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observed cooling as described in the following. The observa t ion shows that the radiative-cool i ng co ronal 
loops are cooling exponentially ( Aschwanden and Terrad^ . 20081 : Morton and Erdelvi . 20091 . 2010i ) and 
the cooling profile has the form 


T = Tojexp(-), 

Tn 


( 6 ) 


where Tg^ is the initial temperature at t = 0 and is the cooling time scale. 

In orde£tomatc^^he_obserwdjxgone^i^^ooling^oft^_b^^gI 2 and£limii.Ei'delyi, Al-Ghafri, and 
Morton (12011 ') assigned thermal conduction to be the essential cause of cooling in their model. Further 
to this, an unsgedfigd^hermod^namksourcefor^creating^he^lggm^^ggling^aUMSSiisii^ Al-Ghafri 


and Erdelvi (|201,‘t l and" Al-Ghafr i et al\ ()2014[) t o investigate longi tudinal MHD Waves in dissipative 
time-dependent plasma. Moreover, Morton. Hood, and Erdelvi ( 2010h assumed that the plasma is cooling 
radiatively and the radiative loss function has the form i5p, where the loss term is assumed to follow the 
Newton cooling. Hence, we assume here that the radiation term p^Q(T) ^ Sp. Therefore the background 
plasma state with no background flow satisfies the equations 


vq = 0, po = const.. 


=4> Vpo = 0, 
Po ^ To, 


Po = —poTo, i 
P 

R dTo _ dTo crj, 

— Po-J— = -OPO, -J— = 
pL at at 


(7) 

( 8 ) 

(9) 


where the 0 index denotes background quantity and <5 is a small quantity. Equation m gives the solution 


To = To* exp{-6t). (10) 

Gomparing Equation (HOD with Equation (|6]) it is obtained that S = 1/tc and this gives a justification 
of takine^Jg^^g^jmal^^hgrgJhg^^^lgrvgijgglingJimgg^^rg^JQQj^^^^^QQQj^Morton, Hood, and 
Erdely i. l2mnh . 

Now, we linearise the governing Equations (di-Eii) about the background state by writing all the 
variables on the form 


f{z,t) = fo{t) + fi{z,t), 

and neglecting nonlinear terms, where the subscript 1 indicates the perturbed quantities. Thus, the 
linearised MHD equations for the longitudinal motion can be reduced to 

d\i jpod\i d^T, di’i 

It is more convenient to use the non-dimensional variables to solve the governing Equation (HU). Hence, 
we introduce the dimensionless quantities 




Vi = 


Vi 


Ti = 


h 

Tf 


Oi 



( 12 ) 


where P is the period of the loop oscillation, Cg is the dimensionless sound speed, Cgi is the initial sound 
speed, and we put P = L/cgi- In what follows we drop the tilde. 

Applying non-dimensionlisation for Equation dill) , we arrive at 


d^vi .^d^vi cr ^ad^Ti dvi 


(13) 
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Here, the dimensionless constants a and e represent the strength of thermal conduction and radiation 
and are defined by 


( 7 - 1 )/ 1 koToY^ 

RLy/'TPoi po 


e = P5. 


(14) 


Both quantities a and e are found to be small under standard coronal conditions in which 7 = 5/3, /t Ri 0.6, 
R = 8.3 X 10^ m^ s“^ and kq ~ 10“^^ m^ s“^ If we take L = 100 Mm and Tqi = 0.6 — 5 MK 

as typical coronal values then we obtain 0.0068 ^ <t < 0.48. 

Now, the governing Equation (USD can be written in the form of wave-like equation as 


d 

[1 

/ d'^vi 


dt 

[cl 

^ dP 

^ dz^ )_ 


ad 9/1 9 ui\ 

7 ^ dt \ Cg dt J 


52 ^ 1 ] d 
dz^ ~'^dt 


Vl dt ) • 


(15) 


Because we are interested in investigating the damping of standing waves so it is necessary to introduce 
the boundary conditions at z = ±1/2. Therefore, as the loop ends embedded in the dense photosphere it 
is suitable to assume that the perturbed velocity vanishes at these ends. 


= 0 at z = ± 1 / 2 . 


(16) 


In case of e = cr = 0 (be. in the absence of radiative cooling and thermal conduction). Equatio n (TTKll 


^qu 

represents a simple wave equation with constant sound speed (see lAl-Ghafri and Erdelvil . 1201311 . The 
effect of cooling and thermal conduction on longitudinal slow MHD waves will be investigated in the next 
sections. 


3. Analytical Solution 

Now, we aim to solve the governing Equation (HID using the WKB method. Thus, we need to introduce 
two slow variables ti = et and a = ed, so Equation HID becomes 


A 

dti 


1 / 2d‘^'vi\ 

dti dz^ J 


a d 


7e 


d ( I dvi 


dti \c1 dti 


ui(z,ti) = ( 3 (z,ti)exp (ie ^ 0 (ti)) , 
where function Q is expanded in power series with respect to e, i.e. 

Q — Qq ± c Qi ±. •. • 

Substituting Equations dUD and (HI into Equation (ED and taking the largest order terms in e (order 
of e“^) we obtain 


d^vi 1 

2 d 

/ 1 dvi\ 

(17) 

9z2 _ 

^ dti 

Ul dtij ■ 

s the form 






(18) 




(19) 


^ + i/<3„ = o, 

where uj = d 0 /dti and Qq is subject to the boundary conditions 


Qo = 0 


at 2 = ±^. 


( 20 ) 


( 21 ) 


according to Equation ED- The general solution to the boundary-value problem. Equations ED and 
ED, can be given by the Fourier series in the form 


Qii{z,ti) = ^ A„(ti)cos((2n± l) 7 rz), = (2n ± l) 7 rCs, n = 0,1,2, - • •. 


( 22 ) 


n—0 
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In this study, we are only interested in the fundamental longitudinal mode corresponding to n = 0. Thus, 
Equation (1^ reduces to 


Qo{z, ti) = A{ti) cos (ttz) , 


(23) 


where A = Aq and uj = ujq. Function A{ti) refers to the amplitude of standing oscillation. In order to 
find the function A{ti) we collect terms of the order of e“^. Then, we obtain the equation 


d'^Qi ^ j_ 

c? ^ c? 


dw 

I_ 

’dti 


3-r)^ + auj^c^ 


Qo + 3a; 


dQo 5 cl d'^Qo cl d^Qo _ a cl d'^Qp 
dti 2 uj dz"^ uj dtidz'^ 7 uj dz^ 


where the function Qi satisfies the boundary conditions 


Qi = 0 


at z = ±- 


(24) 

(25) 


The Sturm-Lioville problem, Equations (1241) and (1251) , has a non-trivial solution when the right-hand side 
of Equation (1241) satisfies the compatibility condition. To obtain the compatibility condition we multiply 
Equation (1^ by Qo, integrate with respect to z over [—1/2,1/2] and use the integration by parts. Hence, 
we arrive at 


r. 1/2 


1-1/2 Cs 


9 dc(j 

—UJ 3-|^—h auj^cl ) Qq 3uj Qo 


dQo , 5 cl ^ d^Qo 
dh 2uj dz^ 


fin 

uj^°dtidz^ 


CTcI d^Qo 
70 ;^° dz* 


The solution to Equation (IMl) . with the aid of Equation 
the form 


dz = 0.| 
(26) 

gives the amplitude of standing wave in 


A{ti) = A{0) exp [clih) - 1]^ . (27) 

This Equation in terms of the scaled variables can be written as 

A{t) = H(0) exp ^ “ 1]) > (28) 

where cl{t) = exp(—5et/2). In case of e = 0 it can be rid of the quantity e ,in the denominator of second 
term in the exponent, using Taylor expansion for the sound speed. 


4. Numerical Evaluations 

In this section, the analytical solution describing the temporal evolution of longitudinal standing-mode 
amplitude are studied using numerical evaluations. Typical coronal conditions are exploited to calculate 
the wave amplitude numerically. Then, the results are depicted to show the behaviour of MHD waves in 
radiatively cooling coronal loops. 

4.1. The effect of radiative cooling 

In the absence of thermal conduction {a = 0) Equation (l2^ reduces to 

A{t) = H(0) exp • (29) 

It is clear that from Equation (1291) the amplitude of oscillating loops decreases with time due to radiative 
cooling. To give more insight into amplitude variations we take e S [0,0.5] as typical values for solar 
corona. Figure [I] shows that the cooling causes a strong damping for the coronal loops. This result is 
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Oscillations Amplitude 



t 


Figure 1. The amplitude of the standing wave with different values of e (0.0, 0.1, 0.3, 0.5) representing the effect of the 
cooling on the amplitude of standing wave. The time is measured in units of L/c^i. 


applicable to TRACE loops with temperature Tq = 1 — 2 MK where radiation is the dominant damping 
mechanism. 

In contrast with the result obtained by lAl-Ghafri and ErdelvH (2o3), who assumed the cooling by 
unspecified thermodynamic source, the cooling by radiation brings about an attenuation in the amplitude 
of waves where a strong cooling (e = 0.5) leads to a strong damping. 


4.2. The eS'ect of thermal conduction 


In the absence of cooling (e = 0) Equation (1^ . after using Taylor expansion for the sound speed, becomes 


A{t) = A(0) exp ( ^ ^ 


(30) 


This expression is consistent with its counterpart in I Al-Ghafri and Erdelvil (l2013f ). The amplitude of 
oscillation is damped and the strength of damping depends on the value of thermal conduction, cr. Taking 
into account that thermal conduction, a, is calculated in an order of e to obtain the amplitude expression, 
so Equation dsni) is physically applicable for a small a. This means that Equation dsni) determines 
the influence of weak thermal conduction, tr <C 1. The variations of amplitude are studied on initial 
temperatures Tq = [0.6,1, 3, 5]® K which correspond to cr = [0.0068, 0.019, 0.17, 0.48]. 

In Eigure[2]we present the effect of varying the magnitude of the thermal conduction, cr, on the damping 
rate of the standing acoustic wave. As we can see that the increase of thermal conduction gives rise to a 
strong decline in the amplitude of standing oscillations due to the presence of thermal conduction. This 
result is mostly expected because thermal conduction is the essential cause of damping for the observed 
hot coronal loops, especially in the region of temperature Tq > 3 MK. 


4.3. Combined effects of radiative cooling and thermal conduction 


Now, we investigate the combined effects of radiation and thermal conduction on damping the amplitude 
of standing slow magneto-acoustic oscillations in radiatively cooling coronal loops using Equation (l28ll . 
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Oscillations Amplitude 



Figure 2. The amplitude of the standing wave with different values of cr (0.0068,0.019,0.17,0.48) representing the effect 
of thermal conduction on the amplitude of standing wave. The time is measured in units of Ljcsi- 


The temporal evolution of longitudinal standing-mode amplitude for various values of e and the initial 
loop temperature is exhibited in Figure [3] We can see a remarkable changes in the wave amplitude for 
different temperature regions. For example, Figures [3^1 and [3 )d indicate that the emergence of cooling 
enhances the rate of damping caused by thermal conduction for loops with temperature Tq < 1 MK. 
On the other hand, the damping rate of wave amplitude starts to decrease due to the plasma cooling 
when the loop temperature approaches 3 MK onwards and the reduction in damping grows quickly by 
increasing the amount of cooling. However, the state of damping is still rising strongly with time in the 
absence of cooling as displayed in Figures (3]: and[3Ji. 

Overall, t he scenario of damping in loops with temperature more than 3 MK is in agreement with that 


Uverall, t ne scenario ot damping : 
obtained bv lAl-Ghafri et al\ ( 20141 ) 


5. Discussion and Conclusion 

In this paper we have investigated the damping of standing longitudinal MHD waves due to radiation 
and thermal conduction in cooling coronal loops. The plasma cooling is assumed because of radiation 
method. The radiative function is postulated to have the form dp to match the observed cooling which has 
approximately an exponential profile. Subsequently, the temperature in the loop decreases exponentially 
with the characteristic time scale which is much longer than the characteristic oscillation period. The 
assumption of the low-beta plasma reduces the MHD equations to one-dimensional system for standing 
acoustic waves. We have used the WKB theory to obtain an analytic solution for the governing MHD 
equations. 

Typical coronal conditions are applied to evaluate the evolution of amplitude with time numerically. 
The results show that the radiative cooling enhances the damping rate of coronal loops with temperature 
To < 1 MK while in the region of temperature Tq > 3 MK the damping is reduced gradually by cooling. 
In comparison with radiation mechanism, it is found that thermal conduction is not sufficient to cause 
a strong damping for very cool loops. However, the damping of hot coronal loops are mainly dominated 
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by thermal conduction, where the amplitude of standing slow MHD waves experiences a rapid damping 
with time in the lack of cooling. 

According to the results, the rate of damping of coronal oscillations seems to increase continuously 
with time until reaching its maximum in the region of temperature 4 < T < 6 MK, and then decreases 
onwards. Eventually, slow magnetosonic waves will propagate almost without damping when thermal 
conduction approaches infinity. 
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